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Abstract. - We propose a simple dynamical model for a glass transition. The dynamics is 
described by a Langevin equation in a piecewise parabolic free energy landscape, modulated by 
a temperature dependent overall curvature. The zero-curvature point marks a transition to a 
phase with broken ergodicity which we identify as the glass transition. Our analysis shows a 
connection between the high and low frequency response of systems approaching this transition. 



The glass transition in both spin systems and supercooled liquids is heralded by anomalously 
slow relaxations ^. Mode-coupling theory provides a framework for analyzing the dynamics 
of the supercooled state Q and has been remarkably successful in predicting the general features 
of relaxations in a highly viscous supercooled liquid. However, detailed comparisons with 
experiments indicate that a complete description of the glass transition is still lacking. An 
alternative description of glass formation relies on the presence of multi-valleyed free energy 
surfaces!^ and the activated dynamics resulting from the presence of these traps [Q. 

Recent experiments indicate that the approach to the glass transition has some universal 
features [H when viewed in terms of the frequency-dependent response of the system. In both 
supercooled liquids j5| and spin glasses ||^, a frequency-independent behavior at high frequencies 
and the presence of three distinct regimes in the relaxation spectrum seem to be generic features 
of the approach to the glass transition. This behavior is different from that observed near a 
critical point where the high-frequency response is thought to be uninteresting [Q. In this paper, 
we present a dynamical model which provides a good description of the frequency-dependent 
response near the glass transition. 

Dynamics of systems approaching the glass transition have been modeled by random walks 
in environments of traps Q]. What distinguishes our model is the presence of an overall 
curvature modulating the landscape of the traps, and the description of the dynamics within 
the traps. The curvature is introduced to model the coupling between a set of variables 
which are inherently frustrated and another which tends to remove the degeneracy leading 
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to the possibility of a unique global minimum and many local minima of the free energy |Q. 
An example of such variables could be the spins and strain fields in a spin glass |^ or a 
frustrated antiferromagnet , or a local orientational order parameter in a liquid coupled to 
local distortions of the liquid [|l0|. The vanishing of the overall curvature is identified in our 
model with the glass transition. This is reminiscent of models where the glass transition is 
associated with an instability [Tl|| or an avoided critical point ||T^ . 

The presence of the curvature (i?) ensures that an equilibrium distribution exists for all 
R > 0. At i? = 0, this distribution becomes non-normalizable and there is a transition to 
a phase with broken ergodicity^, As this glass transition is approached, the relaxation 
within the traps leads to the appearance of a regime in frequency space, R << lu < a/R, 
where the response is essentially frequency independent. The parameter a defines how fast R 
approaches zero as the glass transition temperature is approached. At frequencies smaller than 
a cutoff, ~> as i? "> 0, the response is described by a power law with a positive exponent. 
The exact relationship between ujc and R and the exponent of the power law depend on specific 
characteristics of the model. In particular, we find that if there is no correlation between the 
depth of a trap and its position inside the megavalley, the peak approaches zero algebraically 
and is accompanied by a flattening out of the low-frequency response. This is reminiscent of 
the observations in spin glasses 1^. If, on the other hand, correlations are introduced such that 
the deeper subwells are positioned further away from the center of the megavalley, then the 
shift of the peak follows a Vogel-Fulcher type law and there is no evolution of the low- frequency 
power law: a scenario that is reminiscent of the behavior of supercooled liquids^. 

Fig 1 shows a set of relaxation spectra for the model with no correlations. It demonstrates 
clearly that the low frequency power law approaches zero as the curvature is reduced and that 
there is a high-frequency regime, beyond the peak which grows and becomes flatter as the 
curvature is reduced. Whether or not the high-frequency cutoff to the power law {oj ~ a/R) is 
visible within the same window as the peak (w ~ R) depends on the parameter a. The inset 
shows that there are three distinct frequency regimes and these have a striking resemblance 
to observations in spin glasses]^. 

Each trap, or valley, is represented by a subwell and characterized not only by the time it 
takes to escape from it but also by its internal relaxation time. All subwells in the free-energy 
as well as the overall megavalley are assumed to be parabolic (c/Fig. 2). Each subwell is 
parameterized by its curvature r„, width A„, position of the center </)5^ and position of the 
minimum C„. The set of {C„} is fixed by the requirement that free energy F((f>) is a continuous 
function. To simplify the picture even further, we set all A„ = A which then automatically 
fixes = nA. The curvatures are taken to be independent random variables picked from a 
distribution P(r,n). This defines a free-energy functional; 



where /i„ is the support of the n-th subwell. 

The dynamics is modeled by relaxation in this free-energy surface and is defined by the 
Langevin equation 



where 77 is a Gaussian noise with zero average and variance < r]{t)ri{t') >— TS{t — t'). The 
temperature scale is set by = A^/F and it is useful to introduce the "effective" curvature of 
a subwell; i?„ = i?-f r„. Specific features of the distribution P{r,n) affect the detailed nature 
of the response. The assumption that the curvature of each subwell is uncorrelated with its 
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position, P(r, 11) ~ Pir) is the easiest to implement and is the scenario that we examine in 
detaiL The changes occurring due to correlations are then discussed in the light of these 
calculations. 

A natural candidate for P{r) is an exponential distribution P{r) — e^^"^//3Q. This type 
of distribution has been observed in many spin-glass models]^, where it leads to power law 
distributions of escape times. Such distributions have also been observed in some Ising-like 
systems (without quenched disorder) which demonstrate glassy behavior [^j. It has been shown 
that, in absence of the overall curvature, P = Pq is the temperature where the system falls out 
of equilibrium!^. Our model has two characteristic temperature scales, one set by /Sq and the 
other by R. 

Although Eq.(|) is piecewise linear inside each subwell, the long time evolution of the system 
is a non- uniformly biased random walk in a random environment and there is no exact solution 
to this problem. However, the existence of an equilibrium distribution can be proven rigorously. 
The equilibrium probability distribution is given by exp{—f3F{(j))) as long as the integral of 
this function over remains finite [p^. A simple calculation shows that this integral diverges 
as l/(/3o — 0)VR- In the absence of the overall curvature the 1/\/R factor gets replaced by 
the system size and the probability density diverges as /3 ^ /3oQ- As i? ^ 0, the system falls 
out of equilibrium and there is a change from exponentially decaying correlations to power-law 
correlations in time. This transition can be studied by analyzing the equilibrium correlation 
functions and the response functions associated with these through the fluctuation dissipation 
relation 0. 

A non-zero overall curvature and the requirement that the subwells match at their bound- 
aries, lead to a dependence of the effective barrier heights on the position and curvature of the 
subwells (c/ inset of Fig. 2). The wells become asymmetric and the lower barrier of the n-th 
well is given by 

AF„^^{l~np,,f, (3) 

Zpn 2 

where p„ = R/Rn- This expression is true for npn < 1/2, otherwise AF = 0. Clearly, the 
dynamics in subwells AF > and those with AF = is very different. In fact, it is possible 
to talk about internal dynamics only when the internal relaxation time (rn)~^ is shorter than 
the escape time, which is ~ /3exp(/3A_F'„(r)). This suggests writing the correlation function, 
C(<), as a sum of two components: internal relaxation and barrier crossing contributions from 
the wells with non-zero barriers and a "free relaxation" inside the megavalley coming from the 
subwells where the barriers have vanished due to the overall curvature. The weakest point in 
our analysis is the treatment of the free relaxation. In the absence of an exact solution, we 
have described this as simply a Debye response weighted by the effective number of barrier-less 
wells: 

C{t) = C{tf''P + C{ty'''"' (4) 

n 
n 

In writing these expressions, we have used that fact that an equilibrium distribution exists 
and the probability of finding the system in subwell n is given by the Boltzmann factor, 
exp(— /?F™*"), where F™™ is the minimum value of the free energy inside a subwell. Assuming 
that the local curvature and the position are uncorrelated, and estimating the contribution of 
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the trapped and free relaxation from the "average" weights of barrier and barrier-less wells at 
each position n, leads to: 



np„<l/2 ^ Rrt 

where the first term describes the internal relaxation within subwells and the second term 
describes activated processes. The free part can be written in a similar fashion. 

One of the most interesting aspects of our model is the high-frequency response and the 
origin of this can be understood from an analysis of the limiting model where i? = and 
/3 — > /3oQ. In this limit, the hopping term in (7*'""^ becomes identical to the one analyzed in 
[Q. In frequency domain, the imaginary part of the susceptibility, x ('-^) — ^Ciuj), is known 
to grow as u>^^o~^yi^ at small frequencies and decays at large uj as 1/w. However, there is a 
new feature arising from the internal dynamics which drastically changes the high frequency 
behavior. The internal relaxation part of (7*''''^ is given by 

/■oo -{j3o-P)r 

As /? ^ /3o, the contribution of this term to x{'^) becomes (7r/2 — atan(r*/a;)); an extremely 
slowly decaying function of frequency which has a high-frequency cutoff ~ l/(/3 — /3o)[0. 
The total frequency-dependent response, therefore, behaves as u)'^P°~l^^/P for w ^ and is 
a slowly-decaying function at large frequencies. The free (Debye) part is not present in this 
model with i? = 0. 

The overall curvature changes the effective barriers heights and therefore influences the 
response at both high and low frequencies. To simplify the analysis, we adopt a specific 
relationship between R and /3o, and choose the temperature at which i? = to be at /3o in 
such a way that a(/3o ~ P) = PoPR and take a = 1. The equilibrium relaxation spectra are 
then given by. 

In this expression, the first term arises from the internal relaxation, and is identical to (^, 
except for the cutoff being determined by R. The second term is the "hopping" contribution, 
and can be rewritten, after a change of variables, as 

^^^yrap.hop^ J2^-PRn / . (8) 

P ^ Jo + z IP 

The overall curvature enters through the upper cutoff in the integral, and as i? ^ 0, exhibits 
the slow dynamics discussed in the previous paragraph. There is a tradeoff between the hopping 
dynamics and free relaxation with the hopping contribution decreasing as R increases and more 
and more wells become effectively "free" . The free part is given by, 

nR 



1 prLn 



and leads to a Debye spectrum with peak at lu = R. Numerical results based on this description 
of the response are shown in Fig. 1. The contribution from the deep traps shows one essential 
characteristic of the model; the change in the high-frequency response accompanying the 
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low-frequency changes as i? — )■ 0. The inset demonstrates that the free part is responsible 
for the shift in peak frequencies and the appearance of an intermediate regime. 

There are three different regimes of the response: low frequency power law, associated with 
hopping between different subwells; Debye-like peak coming from barrier-less relaxation and 
high frequency power law decay as a result of superposition of multiple single relaxation time 
processes. The weakest point in this construction is the intermediate, Debye contribution which 
in principle overlaps with both low and high frequency tails. This crossover region we expect 
to be very sensitive to all the approximations made. On the other hand, the low-frequency and 
high-frequency regimes are relatively insensitive to these approximations since they depend on 
processes related to deep subwells. 

This picture is very similar to the experimentally observed response in spin glasses 
with the response flattening out at both high and low frequency ends and the peak moving 
much slower than exponentially. In supercooled liquids, only the high frequency piece flattens 
out, and the peak shifts towards zero frequency according to a Vogel-Fulcher law|^, |^. From 
the viewpoint of our model this difference can be ascribed to a difference in the nature of 
correlations in the two systems. In a structural glass, the crystalline state is the absolute 
global free-energy minimum. The subwells of our model correspond to metastable states and 
the megavalley represents the states accessible in the supercooled phase, with the crystalline 
minimum lying outside this region. This suggests a model where the depth of a subwell and 
its position are correlated and the deeper wells are situated further away from the minimum 
of the megavalley. A correlation in P{r,n), of the form suggested above, leads to an upper 
cutoff in the distribution of barrier heights (c/ Eq. (^). This cutoff leads to a maximum 
escape time t^sc ~ exp(l/i?") where a > is the exponent of a power law describing the 
correlation between the depth and the position of subwells. Beyond this time scale, there is 
only barrier-free motion in our model and, therefore, at frequencies lower than ojc = {tesc)~^, 
we predict that xi'-^) ~ ^ there is no flattening out of the low-frequency response. The 
upper cutoff does not have a large influence on the response arising from the internal dynamics 
of the subwells. The high-frequency cutoff is still approximately given by 1/R as can be easily 
seen from Eq. (^). 

In conclusion, we have demonstrated that the basic features of the relaxation spectra near 
the glass transition can be understood on the basis of a multivalleyed free-energy surface with 
an overall curvature which goes to zero at the glass transition. The spectrum crosses over from 
being pure Debye at large curvatures to one with three distinct regimes with an asymptotic 
high-frequency power law characterized by an exponent approaching zero as the curvature 
goes to zero. The model is also characterized by the distribution of curvatures of the subwells. 
Our analysis suggests that the relaxation spectra of spin glasses and structural glasses can 
be described by the same underlying model with different correlations in the distribution of 
subwells. 

The original motivation for constructing this model came from observations in a frus- 
trated spin model without quenched disorder |8| whose phenomenology is remarkably similar 
to structural glasses. Simulations of this model indicated a free-energy surface with an overall 
curvature and the vanishing of this curvature was accompanied by the appearance of broken 
ergodicity and "aging" . This same model also showed a power-law distribution of trapping 
times in the subwells This spin model can, therefore, be viewed as a microscopic realization 
of the dynamical model presented here. Since it can also be viewed as a model for real glasses, 
further studies of this model should lead to an understanding of the connection between spatial 
inhomogeneities and the anomalous relaxations near a glass transition and clarify the relevance 
of our dynamical model to real glasses. 
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Fig. 1. - Contribution of deep traps to the imaginary part of susceptibility, x {^), for the model with 
no correlations, calculated using Eq. (Q), for different values of R. The peak intensities have been 
matched, ujo is a microscopic frequency scale determined by the parameters of the model. The inset 
shows the full relaxation spectra, including the contribution from the "free" subwells for il = 0.1, 0.05 
and 0.01. It is seen that the free part is responsible for the shift of the peak and the appearance of 
an intermediate regime. The long-dashed lines show the high-frequency power law extending up to 
uj ~ a/R. The parameter a = 1 in the main figure but a = 100 in the inset. 




Fig. 2. - Free energy landscape for two different values of R with a fixed distribution of {r^}. The 
inset shows how the overall curvature modifies the heights of the barriers between the subwells. 



